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METHOD AND APPARATUS FOR ANALYZING A DISTRIBUTION 

Field of the Invention 

The present invention relates to an apparatus for analyzing a distribution. In another 
aspect, the invention relates to a method for signal analysis. 

Background of the Invention 

For the last decade, the speed of modern digital electronics has significantly increased. 
To use SONET as an example, its bit rate quadrupled several times from OC-12 (622 Mbps) to 
OC-192 (10 Gbps). Currently, OC-768 (40 Gbps) is on the horizon. To use Intel's Pentium 
family microprocessors as another example, their clock speed was 60 MHz and 66 MHz in 1993 
for the original Pentium processors, and it is now 1.4 GHz and 1.5 GHz for the Pentium 4 
processors. For such high-speed applications, signal integrity has become an increasingly 
important issue, among which, jitter remains to be a critical performance criteria. 

Summary of the Invention 

In accordance with the present invention, the present invention discloses a method, 
apparatus, and article of manufacture for analyzing a measurable distribution having random 
components and deterministic components. The method includes the steps of collecting data, 
constructing a probability density function based on the data such that the probability density 
function is a convolution of deterministic functions and random functions, constructing a 
probability density function based on a deterministic and random convolution model having 
three or more parameters wherein at least one of the parameters are unknown, and determining 
unknown parameters by using a deconvolution process. 

In one embodiment, at least one random parameter is known and deterministic 
parameters are determined based on one or more known random parameters through the 
deconvolution process. 

In another embodiment, all of the random parameters are known and the deterministic 
parameters are determined based on known random parameters through a deconvolution process. 
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In another embodiment, at least one deterministic parameter is known and random 
parameters are determined based on the one or more deterministic parameters through a 
deconvolution process. 

In another embodiment, all of the deterministic parameters are known and the random 
5 parameters are determined based known deterministic parameters through a deconvolution 
process. 

In another embodiment, none of the parameters are known and the deterministic and 
random parameters are determined through a deconvolution process. 

In yet other embodiments, the distribution is a signal distribution. 

WO The invention may be implemented as a computer process, a computing system or as an 

B 

n article of manufacture such as a computer program product or computer readable media. The 

computer program product may be a computer storage media readable by a computer system and 
encoding a computer program of instructions for executing a computer process. The computer 
program product may also be a propagated signal on a carrier readable by a computing system 
f 15 and encoding a computer program of instructions for executing a computer process, 
fy These and various other features as well as advantages, which characterize the present 

u 

V; invention, will be apparent from a reading of the following detailed description and a review of 
j= j the associated drawings. 



III 



20 Brief Description of the Drawings 

The invention may be more completely understood in consideration of the following 
detailed description of various embodiments of the invention in connection with the 
accompanying drawings, in which: 

FIG. 1 is an exemplary illustration of a representative hardware environment for a signal 
25 analyzing system according an embodiment of the present invention; 

FIG. 2 is a flow diagram illustrating the steps performed by the analysis program 
determining deterministic and random model parameters using a deconvolution process 
according to one embodiment of the present invention; 

FIG. 3 is a flow diagram illustrating the steps performed by the analysis program in 
30 analyzing a distribution according to another embodiment of the present invention; 
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FIG. 4 is a flow diagram illustrating the steps performed by the analysis program in 
analyzing a distribution according to another embodiment of the present invention; 

FIG. 5 is a flow diagram illustrating the steps performed by the analysis program in 
analyzing a distribution according to another embodiment of the present invention; 
5 FIG. 6 is a block diagram representing the steps performed by the analysis program in 

using a deconvolution process; 

FIG. 7 illustrates a determination of the PDF for the deterministic component under a 
noiseless case; 

FIG. 8 illustrates a determination of the PDF for the deterministic component under a 
^10 noisy case; and 

5 FIG. 9 illustrates a distribution that has seven total deterministic and random model 

O 

FU parameters according to one embodiment of this invention. 

M. 

gj While the invention is amenable to various modifications and alternative forms, specifics 

*| thereof have been shown by way of example in the drawings and will be described in detail. It 

y i 

b 15 should be understood, however, that the intention is not to limit the invention to the particular 

Li. 

fTi embodiments described. On the contrary, the intention is to cover all modifications, equivalents, 
H' and alternatives falling within the spirit and scope of the invention. 

6 

^ Detailed Description 

The present invention is believed to be applicable to methods and devices for analyzing a 
20 distribution. The examples below are for jitter analysis. While the present invention is not so 

limited, an appreciation of various aspects of the invention will be gained through a discussion of 
the examples provided below. 

FIG. 1 is an exemplary illustration of a representative hardware environment for a signal 
analyzing system 100 according an embodiment of the present invention. A typical 
25 configuration may include a measurement apparatus 102 that measures the time interval between 
two events (start and stop) through counters. A measurement apparatus is disclosed in United 
States Patent No. 4,908,784, which is hereby incorporated by reference. A typical measurement 
apparatus is the Wavecrest DTS-2075, available from Wavecrest Corporation, Eden Prairie, MN. 
Those skilled in the art will recognize that other systems that enable signal/distribution 
30 analysis that are based on real world measurement (i.e., measurements that are non-ideal or 
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subject to uncertainty) would be applicable. Generally, this would include any product that can 
act as a distribution source. These devices include an oscilloscope, Automated Test Equipment 
(ATE), spectrum analyzer, network analyzer, TIA (time interval analyzer), universal time 
frequency counter, and modulation domain analyzer. Other devices may include a CCD, an x- 
ray camera, a MRI, and an ultrasound. 

The measurement apparatus 102 interfaces to a workstation 104 and operates under the 
control of an analysis program 106 resident on the workstation 104. The analysis program 106 is 
typically implemented through data analysis software. One commercially available analysis 
software is the Wavecrest Virtual Instrument (VI) software, available from Wavecrest 
Corporation, Eden Prairie, MN. Other analysis software includes LAB VIEW, MathCad, 
MATLAB, Mathematica, among others. The workstation 104 comprises a processor 108 and a 
memory including random access memory (RAM), read only memory (ROM), and/or other 
components. The workstation 104 operates under control of an operating system, such as the 
UNIX® or the Microsoft® Windows NT operating system, stored in the memory to present data 
to the user on the output device 110 and to accept and process commands from the user via input 
device 1 12, such as a keyboard or mouse. 

The analysis program 106 of the present invention is preferably implemented using one 
or more computer programs or applications executed by the workstation 104. Those skilled in 
the art will recognize that the functionality of the workstation 104 may be implemented in 
alternate hardware arrangements, including a configuration where the measurement apparatus 
102 includes CPU 118, memory 140, and I/O 138 capable of implementing some or all of the 
steps performed by the analysis program 106. Generally, the operating system and the computer 
programs implementing the present invention are tangibly embodied in a computer-readable 
medium, e.g. one or more data storage devices 114, such as a zip drive, floppy disc drive, hard 
drive, CD-ROM drive, firmware, or tape drive. However, such programs may also reside on a 
remote server, personal computer, or other computer device. Computer readable media can be 
any available media that can be accessed by the computing system. By way of example, and not 
limitation, computer-readable media might comprise computer storage media and 
communication media. 

Computer storage media includes volatile and nonvolatile, removable and non-removable 
media implemented in any method or technology for storage of information such as computer 
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readable instructions, data structures, program modules or other data. Computer storage media 
includes, but is not limited to, RAM, ROM, EPROM, flash memory or other memory 
technology, CD-ROM, digital versatile disks (DVD) or other optical storage, magnetic cassettes, 
magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium 
5 that can be used to store the desired information and that can be accessed by the computing 
system. 

Communication media typically embodies computer-readable instructions, data 
structures, program modules or other data in a modulated data signal such as a carrier wave or 
other transport mechanism and includes any information delivery media. The term "modulated 
^10 data signal" means a signal that has one or more of its characteristics set or changed in such a 

B manner as to encode information in the signal. By way of example, and not limitation, 

P 

fy communication media includes wired media such as a wired network or direct-wired connection, 

and wireless media such as acoustic, RF, infrared, and other wireless media. Combinations of 
v0 any of the above should also be included within the scope of computer-readable media. 

y ! 

£ 1 5 Computer-readable media may also be referred to as computer program product. 
r j The analysis program 106 provides for different measurement/analysis options and 

H; measurement sequences. The analysis program 106 interacts with the measurement apparatus 
Q 102 through the on-board CPU 118. In one embodiment, the measurement apparatus 102 

provides arming/enabling functionality such that the apparatus 102 can measure a signal either 
20 synchronously or asynchronously. The signal is fed to the channel input arming/enabling 
controls 120, 122, 124, and 126 to which event that a measurement is made. 
Counter/interpolators 128, 130, and 132 measure the time interval between the start and stop 
events. Interpolators provide fine time resolution down to 0.8 ps. In response to input controls 
120, 122, 124, and 126, multiplexer 134 controls the counter/interpolators 128, 130, and 132 
25 based on a clock 136 signal. Clock 136 is typically a precise crystal oscillator. 

Those skilled in the art will recognize that the exemplary environment illustrated in FIG. 
1 is not intended to limit the present invention. Indeed, those skilled in the art will recognize that 
other alternative hardware environments may be used without departing from the scope of the 
present invention. 

30 The logical operations of the various embodiments of the present invention are 

implemented (1) as a sequence of computer implemented acts or program modules running on a 
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computing system and/or (2) as interconnected machine logic circuits or circuit modules within 
the computing system. The implementation is a matter of choice dependent on the performance 
requirements of the computing system implementing the invention. Accordingly, the logical 
operations making up the embodiments of the present invention described herein are referred to 

5 variously as operations, structural devices, acts or modules. It will be recognized by one skilled 
in the art that these operations, structural devices, acts and modules may be implemented in 
software, in firmware, in special purpose digital logic, and any combination thereof without 
deviating from the spirit and scope of the present invention. 

FIG. 2 is a flow diagram illustrating the steps performed by the analysis program 106 

0 determining deterministic and random model parameters using a deconvolution process 

according to one embodiment of the present invention. The present invention is directed towards 

lj w 

fy analyzing the deterministic and random components of a distribution. In one embodiment of the 

M. 

g, present invention, the analysis program 106 analyzes the jitter of a signal. In other embodiments, 
% 1 the analysis program can analyze noise associated with other physical quantities including 

e 15 voltage, space, pressure, temperature, velocity, current, frequency, wavelength, power, and 

M< 

pj energy. 

H; Jitter in serial data communication is a difference of data transition times relative to ideal 

□ bit clock active transition times. Jitter in digital systems is the difference between an ideal clock 
period and an actual clock period. As in all signals, jitter has deterministic and random 
20 components. Deterministic jitter (DJ) is bounded in its amplitude and can be measured peak to 
peak. Random jitter (RJ) is unbounded in its amplitude and Gaussian in nature. Since RJ is 
probabilistic, it must be quantified by one sigma of standard deviation estimate. RJ is modeled 
by one or more Gaussian distribution. 

Block 200 represents the analysis program 106 collecting data from the measurement 
25 apparatus 102. The data may be physically based or model based. Block 202 represents the 
analysis program 106 constructing a total probability density function (PDF) based on the data 
such that the PDF function defines a distribution. The distribution is typically a histogram, 
although other distributions, including distributions that are defined by variable bin widths or a 
sliding window can be used. The total PDF is typically a convolution of the PDFs of the 
30 deterministic components and random components. 
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Block 204 represents the analysis program 106 constructing PDFs based on a 
deterministic model and random model having a total of three or more parameters. In some 
embodiments there is only one random model parameter (sigma). In other embodiments there 
can be more than one random model parameter. For a signal analysis application, the 
deterministic model parameters can be time locations and magnitudes of the impulse signals. As 
described herein, the parameters are measured in the relative sense. 

In some embodiments there is only one deterministic time location parameter and two 
random model parameters and the magnitudes of the deterministic model are assumed to be 
uniform. In other embodiments there can be more than one deterministic model parameter. 

The deconvolution method preferably has three or more total parameters. The number of 
parameters is calculated using the least number of parameters for the models. For example, a 
three parameter model can have one deterministic time location parameter, one deterministic 
relative magnitude parameter, and one random sigma parameter. A three parameter model can 
also have one deterministic time location parameter and two random sigma parameters. 

FIG. 9 illustrates a distribution that has seven total deterministic and random model 
parameters according to one embodiment of this invention. The deterministic model is 
composed of three Deltas. The magnitude parameters are d(2) and d(3) where rf(l) = 1. The 
time parameters are t2 and t3 where ti = 0. Thus, the four deterministic parameters are d(2), d(3), 
t 2 and t 3 . The random model is composed of three normalized Gaussian functions, each has a 
standard deviation of a l9 cr 2 and cr 3 . Thus, the three random parameters are a l9 cr 2 and cr 3 . 

Block 206 represents the analysis program 106 determining the parameters for the 
deterministic model and random model by using a deconvolution (or decomposition) process. 

A total distribution is a convolution (or composition) of the deterministic distribution and 
the random distribution. In a jitter example, total jitter (TJ) is a convolution of DJ and the RJ. 
Accordingly, the DJ PDF and the RJ PDF together contribute to the TJ PDF. Once the total PDF 
is constructed, a deconvolution (or decomposition) process means to determine the DJ PDF 
and/or the RJ PDF by solving an inverse problem. 

FIG. 6 is a block diagram representing the steps performed by the analysis program 106 
in using a deconvolution process. 

Block 600 represents the analysis program 106 defining an inverse problem. 
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Block 602 represents the analysis program 106 mathematically formulating an inverse 
problem. 

Block 604 represents the analysis program 106 solving an inverse problem. 

The problem can be formulated as, for example, a set of linear equations, a 
convolution/deconvolution, a polynomial division, a transform domain division, and/or other 
known methods. Approaches to solving the inverse problem include optimizer and non- 
optimizer procedures, and iterative and non-iterative procedures. The initial data for these 
procedures can be "unknown" or "partially-known." "Known" refers to the level of knowledge 
of the model parameters. Generally, it is known that deterministic components are modeled as a 
multiple delta distribution and random components are modeled as one or more Gaussian 
distributions. The knowledge can be described as "unknown" or "partially-known" where the 
deterministic or random parameters are unknown or partially known. In various embodiments 
according to this invention the knowledge of the parameters in solving the inverse problem can 
be described as (1) some or all deterministic parameters are known and the random parameters 
are unknown; (2) deterministic parameters are unknown and some or all random parameters are 
known; and (3) both deterministic and random parameters are unknown. 

In the following, we shall illustrate procedures that can be used to solve the inverse 
problem. The skill of art will recognize that the limited number of exemplary procedures shown 
is not intended to limit the present invention. Indeed those skilled in the art will recognize that 
other procedures can be used without departing from the scope of the present invention. 

Let's denote the DJ PDF model as P DJ {t) , where: 



P D A0 = itd(n)S(t-t n ) (1) 

71=1 

InEq.(l), rf(l) = U,=0. 

Clearly, t n represents the location of the nth Delta impulse of the DJ PDF, and 
P DJ (t„) = d(n) whose associated RJ PDF is denoted as, 
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(t-i,) 2 

P* J -A» = - 7 ±-e (2) 



/, , t N represent the first (minimum) and the last (maximum) locations for the D J PDF, 
respectively. There is total N Delta impulses to describe the DJ PDF with N - 1 independent 
relative magnitude parameters (not including d{\) = 1 ) and N - 1 time parameters (not including 
ti = 0). There are a maximum of N standard deviation parameters to describe the RJ PDF. 

Let's denote the TJ PDF model as P TJ (t) . It is clear that the relationship between TJ PDF 
and DJ and RJ PDFs can be expressed as 

p„(t) = i^(«)^_„ 

= £ d(n) g -ffi (3) 

i.e., the TJ PDF is a weighted and shifted sum of multiple Gaussian functions. The measured TJ 
PDF is denoted as P TJ _ M (t) . 



PROCEDURE 1: 

This is the traditional optimizer procedure as is known in the art, which can be used when 
some or all of the deterministic and/or random model parameters are unknown. 

The following illustrates the typical steps involved in this procedure. 

Step 1: Assume an initial condition, i.e., for all unknown Delta impulse magnitudes and 
all unknown Gaussian function standard deviations, we assume certain initial values. 

Step 2: Evaluate the TJ PDF model based on Eq. (3) using these initial values. 

Step 3: Define an error function based on measured TJ PDF and the calculated TJ PDF, 
we will minimize this error function as our goal by manipulating the model parameters to 
achieve a certain accuracy goal. 

Step 4: The procedure can be iterative or adaptive. Once the goal is achieved, the model 
parameters are determined. 
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For more detail on the optimizer procedure, please refer to references, e.g., Optimization 
Theory and Applications, DA, Pierre, Dover Publications, Inc., New York, 1986 and Numerical 
Recipes in C, WH Press, SA Teukolsky, WT Vetterling and BP Flannery, 2 nd , Cambridge 
Univerisy Press, 1992, both of which are herein incorporated by reference in their entirety. 



PROCEDURE 2; 

This is a recursive non-optimizer based procedure that can be used when some or all 

parameters are unknown. 

Assume that we further divide the interval between adjacent DJ samples into M bins, 
and define 



M 



Step 1: Consider the first Delta impulse of DJ PDF, i.e., t = t r We have 



d{\) 



d{\) 



pr— exp r- 



< At * 



(4) 



= P TJ . u (t l -At) 



which yields 



At 



V2 



ln- 



1/2 



(5) 



In 



P TJ _ M (t x -At) 



1/2 



Step 2: Consider the second Delta impulse of the DJ PDF, i.e., t = t 2 .We have 
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d(l) ( At V d(2) ( NAt ^ 
exp — — +-;= — exp 



■jlna x y-yjla^j 42na 2 {^2 a 

( V 
d(l) NAt 
exp -7=^— 

yJ27ta i ^V 2cr iJ 



= P TJ . M (t l+ At) 



2 J 



(6) 



r • r TJ-M \ l 2J 



7T<J, 



which yields 



NAt 



°2 = 



41 



In 



Pfj-M-PRj-^At) J 



(7) 



d(2) = 42^<T 2 [P TJ . M (t 2 )-P RJ . 1 (NAt)] 



Step 3. Consider now the n th Delta impulse of the DJ PDF, i.e., t = t n .We have 



7VA/ 



4i 



In 



(8) 



din) = J2^a„ [P TJ _ U (t K ) - ( ATA/)] 



If P TJ _ M {i) is in a discrete form, P TJ . M (t a ) can be estimated by linear or polynomial 
interpolation. It is advantageous for the time spacing of the DJ Delta impulses to be much larger 
than the RJ sigma values. The process can be repeated for various sampling combinations 
between two adjacent Deltas and take the average values. The numbers of Deltas and bins do not 
have to be uniformly distributed. 



PROCEDURE 3 

This is a non-optimizer based and non-recursive closed procedure that can be used when 
all of the deterministic parameters or all of the random parameters are known. 
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In other embodiments, the closed procedure can include, but is not limited to, solving 
simultaneous equations. 

FIG. 3 is a flow diagram illustrating the steps performed by the analysis program 106 in 
analyzing a distribution according to another embodiment of the present invention. Block 300 
represents the analysis program 106 collecting data from the measurement apparatus 102. The 
data may be physically based or model based. 

Block 302 represents the analysis program 106 constructing a probability density 
function (PDF) based on the data such that the PDF function defines a distribution. The 
distribution is typically a histogram, although other distributions, including distributions that are 
defined by variable bin widths or a sliding window can be used. The PDF is a convolution of 
deterministic functions and random functions. In this embodiment, the random functions have 
known parameters. The random parameters can be one sigma or multiple sigmas. 

Block 304 represents the analysis program 106 determining deterministic function 
parameters based on the known random model parameters by using a deconvolution process. In 
some embodiments there is one known random parameter (sigma). In other embodiments there 
can be more than one known random parameter. In various embodiments, the deterministic 
function parameters include impulse signal time locations and magnitudes. 

FIG. 4 is a flow diagram illustrating the steps performed by the analysis program 106 in 
analyzing a distribution according to another embodiment of the present invention. Block 402 
represents the analysis program 106 generating a jitter histogram by calculating the local 
statistics in each time bin. A histogram is a statistical representation of the distribution of 
measured physical parameters. The bin size may be fixed or variable for a given distribution. 
Mathematically, a histogram indicates how the number of measurements change over the 
measured parameters. The general procedure for generating a histogram is: 1) measurement 
apparatus 102 measures a signal parameter (for example, period, frequency, duty cycle, etc.) 
repeatedly to obtain a statistical sample; 2) analysis program 106 sorts the data in a descending 
(ascending) order; 3) analysis program 106 automatically defines bin sizes, and the total number 
of measurements that fall into bin ranges is then accumulated; 4) A data set of the number of 
measurements versus the various measured parameters, namely histogram, can be composed and 
plotted graphically. Although the present embodiment describes a histogram, the present 
invention may apply to any kind of distribution. For example, amplitude versus time 
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(waveform), amplitude versus frequency (spectrum), time versus time (jitter time function), time 
versus frequency (jitter spectrum). 

Block 404 represents the analysis program 106 searching and determining the tail parts of 
the histogram by the first and second order derivative method. The tailparts may be found in any 
isolated area of the distribution. In the exemplary embodiment, the tailpart distributions start (or 
end) from the edge of the far left (or right) of the histogram to the first (or last) maximum. 
Further information can be found in United States Patent No. 6,298,315, which is herein 
incorporated by reference in its entirety. 

Block 406 represents the analysis program 106 starting the )? (chi-squared) method for 
fitting the tails of the histogram distributions £ is defined as: 



i-t Ay, 

where j; mod is the model expected value as d efined as: 

y,-y ffi (,0) 

^mod y max 

where y max is the maximum value, // is the mean, and a is the standard deviation of a Gaussian 
function. They are the fitting parameters. Jt, and y i are pairs of data which composite a 
distribution. (In the case of the histogram, x i is the measured parameter, 
y i is the accumulated events corresponding to *.). Ay. is the error of y i data. The best fitting 
parameters will be obtained by minimizing the y mod can be any arbitrary function. £ 
provides accurate model parameter deduction using measurements which are subject to errors 
and statistical fluctuation. 

Block 408 represents the analysis program 106 obtaining the Gaussian distribution 
parameters. The parameters // and a are obtained for the first (left) and second (right) tails of 
the distribution. 

Block 410 represents the analysis program 106 constructing a RJ PDF. 
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Block 412 represents the analysis program 106 obtaining the PDF for the deterministic 
component by solving an inverse problem. In one embodiment, a direct time-domain method is 
used to solve the inverse problem. Those skilled in the art will recognize that other methods can 
be used without departing from the scope of this invention. For example, an optimization 
approach, similar to that discussed above, can be used. 

If P DJ has a finite length of N , and P^ has a finite length of M (M > N), then P TJ 
has a finite length of N + M - 1 . Clearly, N x Unit Interval (UI) can be used to represent the 
peak-to-peak value of the deterministic component. Assuming a single Gaussian RJ PDF, Eq. (3) 
can simply be rewritten as, 

N jV 



Ptj (0 = KD = £ d(n)P u (/ - n) =Z d{n)g(l - n) 



whose matrix form can be written as, 

t = Gd 



(12) 



(13) 



i.e., 



LxN 



Ixl 

'(2) 
' (3) 





0 


0 


0 




g(2) 


2(1) 


0 


0 




g(3) 


£(2) 


SO) ' 


0 














g(N) 


g(N-l) 


g(N-2) ■ 


SO) 


d{2) 


g(N + l) 


g(N) 


g(N-\) ■ 


S(2) 


d(3) 


g(M) 
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g(M-2) ■ 


• £(M-N + 1) 


d(N)_ 


0 


g(M) 


g(M-\) ■ 


• g(M-N + 2) 






0 


0 


g(M) 





(14) 



Here we have selected L to be equal to N + M - 1 , then Eq. (14) represents an over- 
determined problem. Matrix G is a Toeplitz matrix, and it is clear from Eq. (13) that matrix G 
is determined by the PDF of the random component. 

In one embodiment according to the present invention, a linear equation formulation of 
the convolution can be used. A PDF estimation for the deterministic component can be obtained 
by solving Eq. (14). 
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The property of matrix G is important in solving this problem, in particular, whether 
there exists an inverse matrix G' x satisfying, 

G' } G = I (15) 
where / is a unit diagonal matrix, so that we can obtain a true solution to Eq. (13), i.e., 

5 rf = G-'l (16) 

However, it is unrealistic to count on the fact that the inverse matrix of G always exists, 
and most likely, the matrix G is ill conditioned, if not singular. Without considering whether 
matrix G has an inverse or not, a simple pseudoinverse approach can be used. 

Given aLxiV matrix G , let's define a NxL matrix G + that is the pseudoinverse for 
plO G . Assuming the inverse of (G'G) exists, then G + is given as 

IU G + =(CG)~ l G (17) 

Bp where G 1 is the matrix transpose of G . We can obtain, 

d + =G + t (18) 
that is the least-squares solution to Eq. (7), i.e. d + satisfies, 

min||f-Gtf|| (19) 

+; by manipulating d. A pseudoinverse is also called the Moore-Penrose generalized matrix 
M inverse. Eq. (18) is a best-approximated solution to Eq. (13), which is also called the shortest 
length least squares solution, the minimum-norm least-squares solution or the minimum 
Frobenius (Euclidean) norm solution. 
20 If the inverse matrix of G exists, pseudoinverse is simply an expensive way to calculate 

the direct matrix inverse, i.e., in this case, the bottom line is that pseudoinverse can still be used 
to calculate G' x . On the other hand, if the inverse matrix G does not exist, we can always obtain 
G + which has some of, but not all, the properties of G" 1 . 

(G f G) is an autocorrelation matrix and can be ill conditioned. In some embodiments, 
25 white noise can be added to the diagonal elements of the matrix to stabilize the solution, as is 
known to those of ordinary skill in the art. In one embodiment, a constant that is equal to 1% 
(77) of the diagonal value of the autocorrelation matrix (G'G) is added to each of the diagonal 
elements before the matrix inverse. 
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Eq. (18) can be utilized to obtain a stable solution only if there is no noise in the 
measurements. We propose that the raw TJ PDF data should be preprocessed for noise 
reduction before using the above algorithm, for example a sub-space noise reduction as described 
in the following example. For further information, please see references, e.g., "Fully Adaptive 
SVD-based noise removal for robust speech recognition," K Hermus, I Dologlou, P Wambacq 
and DV Compernolle, Proc. EuroSpeech 99, vol. V, Budapest, Hungary. Sept., 1999, pages: 
1951-1954 and Mathematical Methods and Algorithms for Signal Processing, TK Moon and WC 
Stirling, Prentice Hall, Upper Saddle River, NJ 07458, 2000, which is herein incorporated by 
reference in its entirety. 

FIG. 7 illustrates a determination of the PDF for the deterministic component under a 
noiseless case in a jitter example. Plot 700 shows the original DJ PDF and the recovered DJ 
PDF through pseudoinverse, which are indistinguishable. Plot 702 shows the recovered PDF 
after stabilization. Plot 704 shows a single Gaussian PDF for the random component. Plot 706 
shows the TJ PDF. 

FIG. 8 illustrates a determination of the PDF for the deterministic component under a 
noisy case. Plot 800 represents the original PDF for the deterministic component. Plot 802 
shows the recovered DJ PDF after stabilization. Plot 804 shows the noiseless TJ PDF. Plot 806 
shows the TJ PDF with noise reduction. Plot 808 shows the noisy TJ PDF. 

The noise can be generated using a random number generator based on Monte Carlo 
method, and it has a mean of zero and a standard deviation of 10% of the maximum combined 
PDF. In order to use the noise reduction algorithm, the noise variance has to be first estimated. 
In this embodiment, this estimation is performed using the data from the tailfitted segments. 

FIG. 5 is a flow diagram illustrating the steps performed by the analysis program 106 in 
analyzing a distribution according to another embodiment of the present invention. Block 500 
represents the analysis program 106 collecting data from the measurement apparatus 102. The 
data may be physically based or model based. 

Block 502 represents the analysis program 106 constructing a probability density 
function (PDF) based on the data such that the PDF function defines a distribution. The PDF is a 
convolution of deterministic functions and random functions. In this embodiment, the 
deterministic functions have known parameters. In some embodiments, deterministic parameters 
can include impulse signal time locations and magnitudes. 
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Block 504 represents the analysis program 106 determining random function parameters 
based on the known deterministic function parameters by using a deconvolution process. 

It will be appreciated that the figures described above represent a methodology, routine, 
and/or logical flow of program steps which may be implemented to perform the method of the 
present invention. Other programming steps may be included and such specific logical flow is 
intended by way of illustration of a preferred routine. 

Although the invention has been described in with reference to analyzing time jitter, it is 
to be understood that the invention defined in the appended claims is not necessarily limited to 
the jitter examples described herein. As an example, other measurable quantities can be 
analyzed including, but not limited to, voltage, space, pressure, temperature, velocity, current, 
frequency, wavelength, power, energy, and financial quantities. 

The various embodiments described above are provided by way of illustration only and 
should not be construed to limit the invention. Those skilled in the art will readily recognize 
various modifications and changes that may be made to the present invention without following 
the example embodiments and applications illustrated and described herein, and without 
departing from the true spirit and scope of the present invention. 
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